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Abstract 

An absolute measurement of the components of the shear rate tensor 5 in a fluid can be found 
by measuring the photon correlation function of light scattered from particles in the fluid. Previous 
methods of measuring S involve reading the velocity at various points and extrapolating the shear, 
which can be time consuming and is limited in its ability to examine small spatial scale or short 
time events. Previous work in Photon Correlation Spectroscopy has involved only approximate 
solutions, requiring free parameters to be scaled by a known case, or different cases, such as 2-D 
flows, but here we present a treatment that provides quantitative results directly and without 
calibration for full 3-D flow. We demonstrate this treatment experimentally with a cone and plate 
rheometer. 
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I. INTRODUCTION 


Measurement of shear rate in a liquid is of great importance for characterization of fluid 
and surface slip properties. One method for measuring these properties is to examine the 
fluid force or torque on an object (e.g. in a rheometer), but this gives an average indication of 
the slip properties over a large area, and does not directly evaluate the behavior of the fluid. 


Velocimetry (LDV) j^, which can give localized and time-resolved measurements, but all of 
these methods measure fluid velocity. Measuring the shear rate with these methods requires 
measuring a velocity component at a series of closely spaced points, from which the average 
slope is deduced; this requires some measurement time at each location and the region of 
interest must be scanned, so time-resolved measurements of shear are typically not practical 
with these methods. Fnrthermore, since hot-wire anemometry requires holding a probe in 


Direct fluid measurements can be taken nsing hot wire anemometry [l] and Laser Doppler 


m. 


the flow, it may alter the flow and is sensitive to nearby snrfaces 

In this letter we demonstrate the use of Photon Correlation Spectroscopy (PCS) as a 
method of measnring the components of S that does not require any sensor repositioning, 
measures the shear rate directly and with minimal averaging, does not reqnire calibration 
or comparison to a reference, and can be performed to within a few percent within minutes. 


motion jb], and have been used to approximately measure three-dimensional flow in a Conette 


PCS me^asurements have been used to probe turbulent two-dimensional flow j^. Brownian 


Q. 


cell and four-roll mill iTll. The mathematical treatments in prior work have either been 
limiting cases (2-D) or have been approximate derivations with free parameters, requiring 
calibration using a reference. In this letter, we detail our analysis and recent results, which 
yield rapid PCS measurements of the shear rate with remarkable absolute accuracy in 3D 
laminar flow and the beginning of secondary flow in a cone and plate rheometer (CPR) with 
no free parameters. 


II. THEORY 


As with the LDV method, PCS measnrements nse laser light that is scattered from small 
particles that seed the flow |S|. The scattered light is collected and the interference between 
the light scattered from pairs of particles results in the decay of the intensity autocorrelation 
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function. At higher shear rates there is a larger spread of particle velocities and thus a wider 
spread of Doppler-shifted frequencies in the scattered light, which decreases the decay time of 
the autocorrelation function. Unlike LDV, which uses a local oscillator (or multiple beams) 
for a heterodyne measurement, PCS is a homodyne measurement; no additional signal is 
mixed in. 

PCS measures the normalized intensity autocorrelation function G{t) = ~ 1; 

where the angular brackets designate an average over time and / indicates the intensity of 


light scattered from the particles 


B- 


Scattering occurs in a volume containing N scattering 
particles defined by the intersection of the input beam and the field of view of the detector, 
which is a collimating lens feeding into a single mode optical fiber. If the input beam is a 
collimated Gaussian laser beam, then the electric field incident on the particle is 

E^nc^dentAt) = 

= (t)+2CJotg-Uj {t^Zirj {t) 2 ) 


where fj is the displacement from the center of the scattering volume to the particle, 
a = — and b = — are the semi-maior and semi-minor axes of the input beam (allowing for 
the possibility of an elliptical beam) where Wa and Wf, indicate the radius at which the field 
drops to i in that direction, and kg is the input beam wave vector along the axis of the input 
beam (|fco| = We have written this using the tensor Zi = aa^ + blr to simplify and 

generalize the form of the equations, where all vector quantities are column vectors. For a 
circular input beam, Wa = w/, = Wq and 2wq is the conventional intensity ^ beam diameter 
commonly given by laser manufactures. 

When light scatters from a particle illuminated by the input beam, the scattered electric 
field is a spherical wave, and at center of the lens of our detector it can be written as 

EscatteredJ (t) = EinddentJ (t) -^-— - Ei^cidentJ (t) -( 2 ) 

\R-rj{t)\ R 

where R is the displacement from the center of the scattering volume to the center of the 
detector along the optical axis of the detector, ojgj is the frequency of the light scattered off 
the particle and Kj is the scattering vector from the particle (|Kj| = n is the fluid 
index of refraction and A is the laser wavelength). Looking at the equations, we see that 
Kj II (R — rj{t)), so Kj\R — fj(t)\ = Kj • (R — Tj{t)) = Kj ■ R — Kj ■ fj(t). We have also taken 
the approximation that R ^ fj (t) in the final form of the denominator of Eq. El 
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Allowing for the possibility of an elliptical detector, the lens-hber detector system has 
effective major and minor axes c = -^ and d = ^ respectively; Wc and indicate the 
radius at which the electric field of a beam emitted from the fiber through the lens drops 
to ^ in that direction. Only a fraction of the light scattered off each particle is coupled 
in the single mode fiber by the lens of the detector. Thus the field in the center of the 
fiber is Ej{t) oc rjj(t)Escattered,j(t ), where rjj{t) is the coupling efficiency, the fraction of the 
power transferred to the detector by the light scattered from the particle. We find that 
for coupling of a spherical wave to a fiber using a collimating lens, treated with a TEMqo 


Gaussian mode. vEt) oc |e = |e ^(b^ 2 r-i(b| 2 ^ where Z 2 = cc^ + (S' 


r1 


ee 


Appendix® . m The spatial distribution of the coupling was not considered in previous 
work.p, lul Then, the electric field in center of the fiber from a single scattering particle is 

—ii^T R 


Ej(t) oc 


R 


-E. 


incidentj 


{t)e 




( 3 ) 


Combining Eq. |3]and Eq. [U we find that the detected field from a single scattering particle 


IS 






( 4 ) 


where cu' = ujQ-\-u}sj is the Doppler-shifted frequency of the light scattered off the particle, 
kj — ko — Kj is the net scattering vector {\kj\ — 4^sin(0j/2); 9j is the scattering angle for 
the particle), Z — Zi + Z 2 , and Ejj (x Eq ^-^— is the magnitude of the electric field 
which is constant. Because of the form of the equation for G(r), its value is independent of 
Ed . 

We assume that the particles are small enough that the velocity of each particle is the 
same as the velocity of the fluid at that location, that the fluid velocity is time-independent 
during the measurement, and that beam is sufficiently small that the velocity of the fluid 
can be described by a first order expansion about the center of the beam. So the position 
of the particles at some later time t -|- r can be written as 


G (^ + ^) = G d) + VT + Sfj (t)r + rdif,j (r) (5) 

where v is the average velocity of the fluid, rdifji'T) is additional displacement due to diffusion 
in the time r between t and t + r, and S is the shear rate tensor with components Sa/s = 
Then at a time t -|- r Eq. 0] becomes 

Ej(t + r) = Ei:)e“*^n^db+^''r+‘Srj(bT+r*/,j(T)]+iu)'(t+T)g-[rj(t)+j;T]^2[rj(t)+j;T] 
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Returning to the definition of G(r), we have 




E {m)E-(t)Et(t + T)E;(t + t)) 

ijkl 




-1 


( 7 ) 


Terms relating more than two particles do not survive averaging over time if the particle 
distribution is random and the particles are spread over a volume larger than the wavelength 
of light used in all three dimensions jl^. Applying this condition, only terms with pairs 
of particles can survive in the numerator. We also note that the expectation values of the 
different sums are independent. Therefore, 


E {Ei(t)E;(t)E,(t + T)El(t + r)) + E {Edt)Ei(t)Et(t + T)E-(t + r)) 


G(r) — 


ik 


ik 


E (E,m;m 


+ 


E (B™(i)E;(t)> 


E(-EiW£;(i + T)> 


-1 


E{B™WE;(t)> 


-1 


( 8 ) 


For particles inside the scattering volume, there will be very little variation in the scat¬ 
tering vectors, so we can let kj ^ k = — k, where k is the average total scattering vector 

and K is the average scattered wave vector (along the axis of the detector). We see that in 
the sums each term is identical in form and integrated over all space, so we have the sum 
of N identical terms. We also note that displacement due to diffusion is independent of 
the particle position, and so is integrated separately, and is probabilistic with a distribution 
('?■)) = V ^,n e~4D7, where D is the diffusion constant and n is the dimensionality of 

the motion j^, 13]. Terms in the denominator with m ^ n are zero due to the properties of 


Gaussian integrals, and thus the first term in the numerator and the denominator are con¬ 
stant and equal. The second term in the numerator can be easily integrated (see Appendix 
m- Therefore 

G{t) = (9) 


where D = is the diffusion constant, ks is Boltzmann’s constant, T is the temperature 
of the liquid, d is the diameter of the scattering particles, and rj here is the dynamic viscosity 
of the fluid. The constant Gq is included to account for detection efficiency considerations. 
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minimum measurement time and noise; its exact value is unimportant. Note that the mea¬ 
surement is not sensitive to components of S involving the velocity perpendicular to the 
plane containing k, kg and k. 

The average velocity is important only at very high velocities {%FZv ^ {SajsY) Q], where 
particles will be pass into and out of the beam during the course of a measurement. We can 
write the diffusion factor as where q = 2Dk^ in the case of no shear; in the presence 
of shear, q also includes effects such as Taylor Diffusion jn]. Writing the factor in this way 
allows us to easily account for the effect of diffusion while fitting our data without worrying 
about the precise value of q under given conditions. This effective diffusion enhancement is 
visible in all data, explaining the initial non-zero slope of G{t), but does not effect the shear 
measurement. The value of Go is obtained during fitting but is not relevant for determining 
the shear rate. Therefore we may write Eq. |9] as 


G(t) = Goe-'^'e-J 


—qr—\k'^SZ 


( 10 ) 


which can be easily fit once the tensors in the exponent are evaluated for the geometry at 
hand. In contrast to previous theoretical treatment j^, Eq. [IHl can be quantitatively fit 
without the need for scaling or calibration, as the unknown factors {q and Go) are unimpor¬ 
tant and all other factors in the equation can be determined independently. We note that 
this formalism returns the previously reported 2-D equations}^, with corrections for the 
inclusion of the coupling effect, when <S, Z and k are written for that case (ie appropriate 
components are set to 0). 


III. EXPERIMENTAL 

Laminar flow is well understood in a CPR, so it is easy to predict the results of the 
shear rate measurements. We use cylindrical coordinates (p, 0, z) with the origin at the 
intersection of the axis of rotation of the cone and the surface of the fiat plate and ^ parallel 
to the rotation axis of the cone. In laminar flow, if the angle between the cone surface and 
the horizontal is a < 6°, it can be shown, using the continuity equation and no slip boundary 
conditions, that there is only a single non-zero tensor component s = ^ = u/a; u is the 
angular frequency of the cone about the vertical (z) axis[l^. In this limit Eq. [TT)1 reduces to 
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k^s w^w 


b d 2 
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G{t) = Goe-'^^e 


( 11 ) 


Where k^f, = k-<p. Note that G{t) in this case depends only on the size of the beam and the 
receiver in the z direction. For laminar flow we can simply fit Eq. [Ulto find , which 

has no unknown parameters except s. In the event we leave laminar flow, other additive 
terms in the exponent may appear; we do not investigate this behavior in this letter. 

To measure the shear tensor, we direct a Gaussian laser beam from a HeNe laser (6 mW, 
A = 632.8 nm, beam diameter 2w = 0.81 mm) into a cylindrical glass tank and through 
the cone-plate gap of our CPR. The tank is filled with water in which we have suspended a 
low concentration of spheres of diameter d=0.4 fim, ultrasonically dispersed, as scattering 
particles (they are small enough to follow the local velocity; the Stokes number is much less 
than unity (l5|). The concentration of particles is kept low to avoid multiple scattering. 
Typically, cone and plate rheometers are used with liquid only in the gap between the cone 


and plate j^; to minimize errors caused by the excess water outside of the gap, we avoid 
taking measurements close to the edge of the cone. With this incident intensity the photon 
counting rate is of the order of 750 kHz to 1.5 MHz (when shielded from ambient light) at 
a scattering angle of 20.5°. The measurement count rate will increase if unshielded, but as 
this light will be uncorrelated the measurement of shear is insensitive to ambient light levels. 
We drive the cone with a stepper motor, and take measurements by orienting the beam 
and receiver in the horizontal plane, parallel to the plate. The collimating lens receives the 
scattered light, which is coupled into the single-mode optical fiber to an avalanche photodiode 
based single photon counting module (SPCM), so that the direction of the received light is 
well defined. The SPCM signal is fed into an ALV-5000 autocorrelation card that converts 
the intensity I{t) into a (normalized) autocorrelation function G{t). The effective sizes Wc 
and Wd of the detector can be found by measuring the dimensions of a laser beam sent 
backward through the fiber and out through the collimator. 


IV. RESULTS AND DISCUSSION 

An example of PCS data and a fit curve to it is shown in Fig. [2h-. An example of fitting 
shear rate vs cone speed and a comparison of the measured shear rate to the theoretical 
shear rate in the laminar limit is shown in Fig. [2b- Secondary flow begins at a; ~ 2 rad/s 
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FIG. 1. System geometry top view. The incident beam scatters light over a wide angular range 
but only one scattering vector /? is collected by the single mode fiber via the collimating lens. The 
light output is fed into an avalanche photodiode photon counter, then the correlation card. The 
correlation yields the (time-averaged) shear rate s. 



FIG. 2. a) Data and fit for a measurement taken at 0.4 rad/s cone speed at a scattering angle of 
20.5°, with a 2wo = 0.81 mm beam and a 2ws = 1.05 mm detector (on a logarithmic scale vs r^, 
so G(t) is a straight line except near = 0). b) The measured and theoretical shear rate as a 
function of cone speed at a scattering angle of 20.5° for a single run. 

and turbulence begins ai uj ^ 15 rad/s [^. Data sets were gathered by increasing the 
cone speed from 0.2 rad/s to 3.0 rad/s in increments of 0.2 rad/s, with data at each cone 
speed averaged from 20 repeated readings at that speed with an averaging time of 10 s each. 
Five complete data sets were gathered, with a full removal of the sample and realignment 
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between sets, to test the repeatability of the measurement. 

We see remarkably good agreement with the theoretical values of shear rate in the laminar 
regime, where the average slope of the Shear Rate vs Cone Speed graph over the 5 data 
sets was 9.61 ± 0.03 rad~^, a 0.6% error compared to the theoretical slope of 9.55 rad~^. 
Slopes from individual data sets showed a statistical uncertainty of at most ±0.5%, while 
individual slopes varied from the theoretical value by a maximum of ±2%, larger than could 
be accounted for by statistical error alone, indicating additional sources of error in our 
experiment. 

Tests of repeatability indicate that the additional error results from limited repeatability 
in assembly and alignment, including the setting of the scattering angle, measurement of the 
widths of the input beam and the lens/hber detector, determination of the direction of the 
velocity at the measurement point (which is needed to properly perform the scalar product 
• 0), errors vertically aligning the beam and the receiver with each other and the gap, and 
error in the distance of the cone from the plate. Despite these issues, our measurements 
were still highly accurate, indicating that the method is quite robust. 


V. CONCLUSION 

We have demonstrated the efficacy of Photon Correlation Spectroscopy for measuring 
the shear rate in a three dimensional system, with a short individual measurement time. 
We show that this measurement is quantitatively precise, removing the need to scale to 
a reference found in some previous work, and tolerant of minor systematic errors. By 
measuring the shear rate directly and rapidly, we open up the possibility of looking at 
phenomena that cannot be examined utilizing slower methods or methods that depend on 
averaging over large areas. 
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Appendix A: Coupling Efficiency 


In order to determine the coupling efficiency rjj (t) , we rewrite Eq. [2] for any observation 
point. 

^-iK\h^g{t)\ 


^scatteredJ (^5 ^ineidentj (^) ' 


\h + git)\ 


(Al) 


where K is the magnitude of the scattering vector parallel to h + g{t) {\K\ = ^), g{t) = 
R — fj{t) is the displacement from the particle to the center of the collimating lens, and h is 
the vector to the observation point from the center of the lens. We neglect the phase factor 
as it has no spatial dependence. Although the vector K changes direction as a function 
of h, the magnitude is constant. 

The field in a plane parallel to the plane of the lens may be represented by the collimated 
Gaussian TEMqo mode 

Effi) = (A2) 


where Efi is the field magnitude. 

The coupling efficiency is given by the overlap integral 





f E seattered,j{h 1 Ejt(^h^dA. 

A 

2 

I ^scattered,j(^d)dA 

A 

jEj{h)dA 

A 


(A3) 


which indicates the ratio of power coupling between the modes through a surface A to the 
total power in each mode j^]. Because G{t) depends on the ratio of intensities, in all cases 
the constant normalization terms in the denominator are unimportant. 

Turning to the coupling term, for the plane containing the lens we have 

2 


Vjd) 


oc 


/ Escattered,j(h) t)Elt{h)dA 

2 

r ^—iK\h-\-g{t)\ ^ 

f ^ dA 

J ^ 

A 


J \h + g{t)\ 


(A4) 


which is non-analytic. However, since the integral is over this infinite plane, the integral 
must be the same on any parallel surface (so long as the reversed detector beam would 
still be approximately collimated when it passes through that surface), and we are free to 
choose the most convenient one. Thus we evaluate that integral on the surface containing 
the scattering particle, where Escattered,j{h,t) = Eincident,jif)^{h - g{t)), which gives 
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Vj{t) 


oc 




(A5) 


We can also approximately evaluate the integral on the plane containing the lens. Looking 


at the spherical wave, \h + g{t)\ = \Jh? -\- g{ty + 2h ■ g{t) ~ g{t) + 


I h-g{t) 

"igit) 9{t) 


, where we 


have Taylor expanded to 1^* order with g{t) > h; we cannot a priori say that h?/2 <C 
h ■ g{t) (as g{t) and h are nearly perpendicular). If we retain both terms in this expansion, 

Tp (h +\ W e V / „„ 

^scattered,jy ^^ ^ineidentj\^) g(^'^ ) 




OC 


^-iKg(t) 


gif) 

Q-iKg{t) 

1 

gf) 



g 4:g{t)‘^—2ig{t)Kw‘^ g — ‘^i 9 ii)Kw‘^ 


(A6) 


If we do not retain h‘^/2, then we have performed the plane wave approximation and 


-iKi g(t) + 


h-g(t) 

9(t) 


Escattered,jih,t) Eincident,j^)- -- , which givCS 


Vj{t) OC 


Q-iKg{t) 




gif) 

^-iKg{t) 


9{t) 


-TiWcWde 


Ai-9(t)-c) wi+{-At)-d) 


(A7) 


A comparison of Eq. IA51 Eq. IA6I and the numerical integral of Eq. IA4I shows excellent 
agreement when g is small compared to Kwc and Kwd, which is precisely the limit where 
diffraction effects are not large. By contrast, Eq IA7I is very different, showing that the plane 
wave approximation is a poor choice in this case (See Eig. 13)). 

Einally, using the exact solution lEq. lABjl . remembering that g{t) = R—fj{t), and noticing 
that R A c and .R T d by definition, we see that 


Vjd) 0^ 
where Z 2 = cc^ + dd^. 




„-{rj{t)-c)'^ -{rj{t)-d)‘^ 




(A8) 
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FIG. 3. Comparison of the coupling efficiency in water at ^ = 200 mm and A = 633 nm, as 
calculated numerically on the plane of the lens, analytically evaluated on the plane containing 
the particle, and approximately integrated on the plane of the lens. We also show the result 
of evaluating with a plane wave approximation (neglecting h‘^/2) to demonstrate that this is an 
overapproximation. 


Appendix B: Evaluation Of G(t) 


The normalized correlation function is 


G(t) 


(I(t)I(t + T)) 


E {Ei(t)E‘{t)Et(t + T)E;(t + r)) 

ijkl 

2 


E {E,ME‘{t)) 


(Bl) 


Terms relating more than two particles do not survive averaging over time if the particle 
distribution is random and the particles are spread over a volume larger than the wavelength 
of light used in all dimensions [l2|. Applying this condition, only terms with pairs of particles 
can survive in the numerator. We also note that the expectation values of the different sums 
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are independent. Therefore, 


G(r) 


E {Edmmm + T)Ei(t + t))+Y. {Edt)EiH)Edt + r)B-(i + r)) 

ik ik 


E {E,„mm) 

mn 


E (Ej(t)E;(t)) 

3 

2 

+ 

2 

E(-Bj(t)q(< + d> 

3 


Y.{Ejt)Em 

mn 

2 


1 


(B2) 


For the denominator, 


K-fW)t = 


Y,(Er.(t)E'n(t)) 

mn 


Y,{E,u{t)Erjt)) 

m^n 


+ 


Y,{EAt)E-(t)) 


ElT.' 


-2f^{t)Zrj(t) 


e 3 


ElN / 


EiN 


TT" 


det(2Z) 


= Ej)N^ 


TT 


det(2Z) 


(B3) 


where eac 
generality 


1 of the N expectation values is identical and the integral is n-Dimensional for 




The first term in the numerator of Eq. IB2I is identical to Eq. IB3[ Eor particles inside the 
scattering volume, there will be very little variation in the scattering vectors, so we can let 
kj = ko — Kj ^ ko — K = k, where k is the average total scattering vector {k is the average 
scattered wave vector which is along the axis of the detector) and cn' ~ oj', where u)' is 
the average total phase shift. We see that each term in the sums is identical in form and 
integrated over all space. We also note that displacement due to diffusion is independent of 


^ (Ex = Y This can be verified by applying the unitary transform a^a = 1, where this 

diagonalizes A, and then performing the integral in eigenspace. Such a transformation does not change 
the determinant. 
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the particle position, and so is integrated separately. Therefore 

2 

J2{EAt)E;{t + T)) = 

3 

Zvt'^^- iw'r^ik'^Srj{t)T^-2ff {t)Zrj{t)^-2JFZrj{t)T\^ |2 _ 


E%e 


X 


OO OO 

^jpirdifj (r) mm) Zr^ity 

^ —no —no 


(B4) 


where P{fdif{T)) is the n-D probability distribution of the displacement due to diffusion. 
The integrals are evaluated over all possible particle positions and diffusion displacements. 
For the integral over particle positions in Eq. IR41 we have 


^iP'Sfj (t)r g-2rT {t)Zrj (t) ^-2v'^Zrj _ 


_ ^—^k^S(2Z) ^S'^ky g— ^ik’^SvyZZy 


(B5) 


where D = 

SttijcL 


det{2Z) 

The n-D probability density for diffusion is P(r(T)) = '' 

is the diffusion constant; ks is Boltzmann’s constant, T is the temperature of the liquid, 
d is the diameter of the scattering particles, and 7] here is the dynamic viscosity of the 
fluid jb, 13]. So for the integral over the diffusion displacements in Eq. IB4I we have 

OO 




-Dk'^ 


(B6) 


V4:TTDT 

— OO 

Putting these into Eq. EH we have (noting that the N terms in the summation are 


identical) 


j2{E,m;(t+T)) 


= EiN^ _ ^ _ -2Uk-T -V^ZVT- -fk^az-^a^kr- 

^ det{2Z) ^ ’ 


TV" 


—2Dk‘^T^—v^ZvT‘^^—\k'^SZ 


2 |Oo^ g-.- .4.-+B 33^..^ ^ ^ 


(zl )_^ This can be verified by applying the unitary transform a = 1, 
where this diagonalizes and then performing the integral in eigenspace. Such a transformation does 


not change the determinant. 
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Therefore, 


G(r) 


(I(t)I{t + T)) _ 

—2Dk‘^T —v^ZvT‘^ —jk'^SZ ^S'^kr^ 

mw 


(B8) 
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